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ABSTRACT 

Two important avenues into understanding the formation and evolution of galaxies are the Kennicutt-Schmidt 
(KS) and Elmegreen-Silk (ES) laws. These relations connect the surface densities of gas and star formation 
(Sgas and E*, respectively) in a galaxy. To elucidate the KS and ES laws for disks where Sgas^ 10"* M0 pc^^, 
we compute 132 Eddington-limited star-forming disk models with radii spanning tens to hundreds of parsecs. 
The theoretically expected slopes (w 1 for the KS law and w 0.5 for the ES relation) are relatively robust to 
spatial averaging over the disks. However, the star formation laws exhibit a strong dependence on opacity that 
separates the models by the dust-to-gas ratio that may lead to the appearance of a erroneously large slope. The 
total infrared luminosity (Xtir) and multiple carbon monoxide (CO) line intensities were computed for each 
model. While Ltir can yield an estimate of the average Eh, that is correct to within a factor of 2, the velocity- 
integrated CO line intensity is a poor proxy for the average Egas for these warm and dense disks, making the 
CO conversion factor (aco) all but useless. Thus, observationally derived KS and ES laws at these values of 
Egas that uses any transition of CO will provide a poor measurement of the underlying star formation relation. 
Studies of the star formation laws of Eddington-limited disks will require a high- J transition of a high density 
molecular tracer, as well as a sample of galaxies with known metallicity estimates. 
Subject headings: galaxies: evolution — galaxies: ISM — galaxies: starburst — stars: formation 



1. INTRODUCTION 

Galaxies can experience a wide variety of star formation 
phenomena, from the relative calm of star-forming regions 
beaded along a spiral arm, to the fury of a nuclear starburst 
burning through a trillion Suns. Yet, remarkably, the rate 
of star formation in such wildly disparate environments can 
be simply related to the density of gas in the star-forming 
region, E» cx E^g, where E, and Egas, are the star for- 
mation rate (SFR) and total gas (i.e., atomic plus molecu- 
lar) surface densities. This relation, known as the Kennicutt- 
Schmidt (KS) law, has been observed to hold over near ly four 
decade s i n Ep^^s with N sa 1.4 (e.g., Schimdt 1959; Ke nnicuttI 
Il998allbt lYao et all I200I iKennicutt & Evansi i2012). and is 
ultimately related to the efficiency of star formation and 
its connection to both loca l and global time s cales (e.g., 
ElmegreenI fT997| IsilB fT997t IKennicutt) Il998bl; lElmegreesI 
2002; Krumholz & McKee 2005; Krumholz et al. 200^ The 
super-linear slope indicates that the star formation efficiency 
(i.e., the ratio of stellar mass formed to the gas mass in the 
star-forming region) increases at high densities. Measuring 
and understanding the KS law over a wide range of star- 
forming environments is crucial to elucidating comprehen- 
sive theories of star and g alaxy formation (see discussion by 
IKennicutt & Evansll2012h . 

Most of the star formation in the universe occurred long 
ago at z > 1 when the majority of massive galax- 
ies were being assembled (e.g., iHopkins & BeacomI 120061 ; 
iPerez-Gonzalez et al.l 12008 1). With the advent of new sensi- 
tive milUmeter arrays and detectors, the last decade has seen a 
significant advance in the study of the KS law at high re dshifts 
(e.g. Daddi et al. 2010b; Genzel et al. 2010; Ivis on et al. 
120111: iFreu ndlich et al. 2013; see Solomon & V anden Bouti 
(|2005|) and Carilli & Walter (2013) for reviews). Interest- 
ingly, these galaxies, which have both much larger star for- 
mation rates and gas densities than local objects, seem to fol- 
low a similar KS slope, but are vertically offset from the lo- 



cal relation (iDaddi et al.ll2010bl: iGenzel etaP 120101) perhaps 
implying a fundamental difference in the stellar initial mass 
function (IMF) in the two regimes. However, many of these 
results rely on the assumption that the troublesome conver- 
sion factor aco = 5]gas/-^co needed to convert from the ob- 
served CO velocity-integrated intensity, Ico, to Egas is bi- 
modal with one value (aco ~ 4) for galaxies with Milky 
Way-like SFRs, and another value (aco ~ 0.8) for galax- 
ies with much larger SFRs (i.e., high-z galaxies). When 
[Narayanan et al.l (12012 1) recently recomputed the KS law with 
a value of aco that varied continuously with /co, the relation 
was no longer bifurcated between the high and low-z galaxies 
and could be fit over a wide range in Egas with N Ri 2 (see 

also lOstriker & Shettyll201 ll) . 

Alternatively, it has been argued (plmegreenlll997t ISilli 



[T997HDaddi et al.ll2010bl;rGenzel et alj|2010h that a more fun- 
damental star formation law is E* cx (Egas^^)", where J7 
is the orbital angular frequency in the star-forming region, 
usually estimated at the observed outer radius of a galaxy. 
Both high and low-z galaxies seem to follow this Elmegreen- 
Silk (ES) relation w ith n fti 1, independent of the exact as- 
sumptio ns on aco d Daddi et ajj l2010bl; iGenzel et al.l 120101 ; 
although iNaravanan et al. (20121) found that the scatter in the 
ES relation was reduced by using the continuously varying de- 
scription of aco)- In this formulation, a larger star-forming 
efficiency is found in rapidly star-forming galaxies because 
of occurring in more compact regions with shorter dynamical 
times. 

Analysis and interpretation of the star formation laws at 
high redshift are hampered by gaps in the theory of star for- 
mation in galaxies at these redshifts. The interstellar medium 
of these galaxies are expect ed to be predominate l y molecular 
(e.g., IS olomon et al. 19971 iTacconi et all 120 lot iDaddi et all 
l2010at iMagdis et al.ii2012l) have gas f ractions of several tens 



of percent (e.g iGeach et all 120111: ISwinbank et"aLl 120121; 



iNaravanan et al.ll2012bh with significant turbulence and largi 
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pressures (e.g-jDownes & SolomoDll998l:lGenzel et al.l2008l: 
IWisnioskietalJ2011h . Simulating this environment from first 
principles is computationally prohibitive, so analytical mod- 
els have been developed that encompass many aspects of the 
physics in a we ll-defined region of parameter space (e.g., 
[Thompson et al.i 2005: Ostriker& Shetty 2011). In addition, 
as eluded to above, there may be severe uncertainties in re- 
lating the observed quantities (usually, infrared (IR) luminos- 
ity an d CO intensity) to the physical p arameters an d S^a, 
(e.g., LS olomon & V anden Boull l2005"t iFeldmann et all 120121: 
iNaravanan et al. 2012|k Insofar as the IR luminosity cap- 
tures the bolometric output of embedded star-forming re- 
gions, there are well known calibrations between the observed 
IR lu r ninosity and SFR ( e.g., Kennicutt 1998b; Riek e et alJ 
l2009t ICalzetti et al.ll2010i : Murphv et al.ll201 ID . More prob- 
lematic is how to convert from the line intensity of a molecule 
that comprises only a tiny fraction of the interstellar gas to 
an estimate of the entire gas mass. Although the CO lines 
are most commonly used due to its high abundance, it has 
been noted that because of their relatively low critical den- 
sities and high optical depth, t hey are likely very poo r trac- 
ers of high dens ity gas (e.g., Krumholz & Thompson 2007t 
INaravanan et al.ii2008) . Indeed, KS laws derived from a high 
density tracer such as HCN have produced a simple linear re- 
lation between and E^as for both high- and low-z galax- 
ies (e.g., Gao & Solomon 2004a,b), although the conversion 
factor ancN is just as uncertain as the one for CO (see the 
discussion by lGarcfa-Burillo et al.ll2012h . Taken together, the 
combined uncertainties in both the existing theoretical frame- 
work and the observational conversions severely restrict the 
use of the KS or ES laws. 

The most extreme star-forming regions are those with 
Sgas ^ 10"* M0 pc^^. In this situation the opacity of 
the dense and dusty gas is so large that the ISM becomes 
optically thick in the IR, radiation pressure dominates over 
the turbulent pressure, and the disk is said to be Eddington- 
limited. Expected to occur in the cores (i.e., central few hun- 
dred pes) of rapidly star-forming disks, these intense star- 
forming regions are difficult to investigate observationally, al- 
though evidence for such high values of Egas have been re- 
cently in ferred at the centers of some z < 1 post-starburst 
galaxies (Diamond-Stanic et al. 120121) and star-forming sys- 
tems (Garcia-Burillo et al.l l2012 ). Their location also implies 
that these Eddington-limited star-forming disks may act as 
a 'bridge' between the black hole environment and the host 
galaxy, ferrying fuel through the disk toward the central black 
hole and p otentially driv ing nuclear activity C Thompson et al.l 
1^05: Ball antynd 120081) . Thus, predictions of the observa- 
tional signatures of Eddington-limited star-forming regions 
are needed to guide future observational studies of the cen- 
ters of both active and inactive galaxies at all redshifts. 
lOstriker & Shettvl (1201 lb analytically investigated star forma- 
tion in environments with 10^ < Egas ^ 10 M© pc~^ when 
vertical pressure support is provided by supernova driven tur- 
bulence, ^nd_gredictedjharthe KS law should have N — 2 
(see also Thomp son et al.l2005-) . in decent agreement with the 
high-z observations if one assumes that aco is inversely cor- 
related with Iqo- In the Eddington-li mited case the local KS 
law is expected to flatte n to = 1 (iThompson et al.l 120051: 
lOstriker & Shettyll201 Ih . but there have been no predictions 
for the ES law nor work on how clearly these laws will trans- 
late for an observation that encompasses an entire > 100 pc 
star-forming disk and contains a significantly variable SFR 



and Egas. In this paper we use the IThompson et al.l (120051) 
model of an Eddington-limited disk to self-consistently cal- 
culate the radial structure of over 100 nuclear star-forming 
disks, and, through the use of CO radiative transfer calcu- 
lations, we investigate the relationship between observable 
(Lifj and Ico) and physical (E, and Egas) variables in these 
most extreme star-forming environments. 

The next section describes the Eddington-limited star- 
forming disk model, the calculation of our model database, 
and the CO radiative transfer method. Both the KS and ES 
star formation laws predicted by these models are presented 
in Section |3] The observational star formation law (the IR- 
CO relationship) that is predicted by these models and its re- 
lationship to the physical KS and ES laws is described in Sec- 
tion|4] A discussion and concluding remarks are presented in 
Section|5] Throughout this paper the total infra red luminosity 
(Ltir ) is defined as the 3-1100 /im luminosity (iMurphv et al.l 
1201 ll) . and helium is not accounted for in Egas or aco ■ Unless 
otherwise specified all CO intensities or luminosities refer to 
the J = 1 — rotational transition. 

2. CALCULATIONS 
2.1. The Eddington-Limited Star-Forming Disk Model 

The Eddington- hmited starburst disk m odel utilized here 
was developed by IThompson et al.l (120051) . and a short sum- 
mary of the relevant details is presented below (see also 
Thompson 2009). 

Star formation is modeled simply as an energy source at 
a location r within a one-dimensional single phase medium 
that rotates with angular frequency in the potential of 
a galactic bulge (modeled as an isothermal sphere) with 
dispersion a and central black hole of mass A/bh (these 
two quantities are assumed to be related by Mbh = 2 x 
10*(o-/200 km s-i)-^ M©; e.g., Tremaine et aDI^OOl). Gas, 
radiation and turbulent pressure driven by supernova all com- 
bine to provide the vertical support against gravity: 

Pgas + eE*c(^+e) =p/i'f^', (1) 

where pgas — pkT/mp, and ph^il^ is the total pressure re- 
quired for hydrostatic balance in the bulge potential (with h 
denoting the gas scale height and p indicating the gas den- 
sity of the star- forming disk; the self-gravity of the disk is 
neglected here). IThompson et al.l ( l2005i) show that in the op- 
tically thick limit the radiation pressure is eE*CTy/2 where 
e is the IMF-dependent efficiency between SFR and lumi- 
nosity (i.e., L — eM^c^, where ill* is the local SFR) and 
TV — KEgas/2 is the vertical optical depth in the IR (k(T, p) 
is the Rosseland mean opacity of the dusty gas at tempera- 
ture T and density p). The turbulent pressure driven by su- 
pernova can also be written as proportional to the local SFR 
density, e^S^c, where ^ « 1 (Thompson et al. 2005) . From 
above, a rms speed Cs can be defined using the total pressure 
as Cs = hft, where it is understood that gas, radiation and tur- 
bulent pressure all contribute to this rms speed. As described 
below, the models studied in this paper all have a radially- 
averaged Ty > 1 and are thus dominated by radiation pres- 
sure for almost all radii. 

The star-forming disk is assumed to be always just unstable 
to gravitational instabilities; that is, Toomre's Q parameter is 
set equal to one. Material with gas fraction /gas is fed onto to 
the disk at a radius rout from the black hole, and, through a 
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hypothesized global torque (pr ovided by, e.g., a spiral insta- 
bility or a b ar; Goodman 20Q3: lEnglmaier & Shlosmanll2004l 
iMacieiewskii i2004) . can slowly accrete towards the center 
with a radial velocity equal to a fixed fraction m of the lo- 
cal rms speed (iGoodm an 2003j), i.e., M = 27rrEgas»TiCs. At 
any radius r. 'Thomp son et al.i(i2005h showed that M and fgas 
can be related through 



/g 



23/2 AfG 



1/2 



(2) 



As the gas moves through the disk, it is also producing stars at 
a rate of E*, which, in the optically thick limit, is determined 
by requiring radiation pressure in the IR to support the disk 
(eq.[T]l as well as Q 1: 



\/2/gasQcr' 



£KCr 



(3) 



where r is assumed to be large enough that the black hole 
potential is negligible. The star formation slowly consumes 
the gas as r decreases, reducing both M and /gas, i.e.. 



M = Af o 



(4) 



where A/out is the mass accretion onto the disk at Tout- How- 
ever, the (5 = 1 assumption causes a rising density in the inner 
regions of the disk that substantially increases Ty so that even 
a moderately increasing can still support the disk through 
radiation pressure. The solution at radius r is determined by 
finding the minimum central gas temperature T and associ- 
ated K to ensure vertical pressure balance, either through the 
sum of gas and turbulent pressure, or, if there is enough star 
formation, radiation pressure. Eventually, either the gas avail- 
able for star formation becomes too small to maintain Q = 1, 
or energy release through gas accretion pushes Q > 1, and 
the calculation ceases. 

Substituting eq. |2]into eq. |3] results in predictions for both 
the ES and KS laws in the disks: 

^* = ' — r^Sgasf^ (5) 



and, using Q to eliminate Vl (see [Thompson et aI]|2005L eq 

3), 



■^gas- 



(6) 



These equations have to be viewed carefully because, as seen 
in eq.m Egas is not an independent variable; its value at some 
radius r will depend on E, at larger radii. As seen in Sect. |3] 
this property may affect the predicted slopes of the star forma- 
tion laws. In sum, radiation pressure supported star-forming 
disks predict a KS relation with N k, 1, and an ES law with 
n « 0.5. Of course, both of these predictions rely on the as- 
sumptions that (5 = 1 and the disk is optically thick in the IR 
and therefore radiation pressure supported. 

In order to determine how these laws will manifest in ob- 
servational surveys of several galaxies, the full disk mod- 
els need to be calculated over a wide range of parameters. 
There are five input parameters for each model: the black 
hole mass A/bh (which determines a), the angular momen- 
tum parameter m, the outer radius rout, the gas fraction at the 



outer radius /gas, and a dust-to-gas multiplicative factor for 
the opacity to account for the e nhanced metalhcity observe d 
in the centers of galaxies (e.g., iMu noz-Mateos et al."2009l). 
As in Thompson et al. (20051), the[Semenov et al. (2003) cal- 
culation of the Rosseland mean opacity for dusty interstel- 
lar gas is used for all models. Calculations are performed 
for each permutation of log(AfBH/Mo) = 7,7.5,8,8.5, 
m = 0.0075,0.01,0.025,0.05,0.075,0.1,0.2, rout = 
50, 100, 150, 200, 250 PC /gas = 0.1, 0.5, 0.9 and dust-to-gas 
ratios equal to 1 x, 5x and 10 x the local ISM value, resulting 
in a suite of 1260 models. While the models are limited to 
font < 250 pc, the results on the star-formation laws do not 
depend on this size (as is seen for the observed KS and ES 
laws), and are valid for any Eddington-limited star-forming 
region. The range of inpu t gas fraction s is consisten t with 
the ob served estimates (e.g.. lTacconi et al . 2010; Magdi s~et al.l 
120 12h that show an increase to ^50% at z ~ 2 (the /gas input 
to the model drops steadily through the disk as the star for- 
mation uses up the available gas; eqs.|2]&|4]i. We emp loy a a 
standa rd Kroupa/Salpeter IMF with ajMF = —2.35 (S alpeterl 
1955; iKroupa & Weidneri l2003) between 1 and 100 Mq and 
aiMF = — 1.3 between 0.1 and 1 M© (where the number of 
stars with masses between M and M + dM is proportional 
to Af "■">'). This IMF is then input into a Starburst99 stellar 
synthesis model dLeitherer et al.lll999) to compute the bolo- 
metric luminosity of a stellar population with a constant SFR 
and with the given IMF as a function of time. The value of 
e = LjMf? is then measured at an age of 10^ yr yielding 
e = 7.1 X 10""*^ (very si milar to value of 6.2 x 10^"* used by 

lOstriker & Shettvll201l]Pl . 

As described in detail by [Thompson et al.l (120051) . roughly 
two different types of star-forming disks result from this the- 
oretical description. In about one-third of the models the gas 
maintains a large optical depth over the entire disk, resulting 
in the temperature in the inner parsec of the disk exceeding the 
dust sublimation temperature and causing a significant burst 
of star formation in order to maintain (5 = 1. These parsec- 
scale bursts could puff up the disk, and, if they do not use up 
all the gas, may both fuel and obscure the central black hole. 
The observati onal pr operties o f these disks have b een stud- 
ied elsewhere (iBallan tyne 2008: lArmour & Ballanty ne 2012), 
and, as they are not directly related to the study of general star- 
forming galaxies, the description of their star formation laws 
are deferred to later paper The remaining models reach the 
maximum SFR, SFRmax, at rout with the SFR declining to- 
ward smaller r. Many of these disks are actually gas-pressure 
dominated or do not extend very far before E* becomes too 
small to maintain (5 = 1. As our interest here is studying ra- 
diation pressure dominated star-forming disks that can extend 
tens to hundreds of parsecs, we selected only models with a 
radially-averaged ry > 1 and an inner radius < 0.04rout for 
further investigation (including less extended disks only in- 

' The assumption of Q = 1 and the location of the star formation 
in the galactic nucleus results in significant gas densities with p rising 
from ~ 10-20 g cm"^ (148 Mq pc' ^ ) to > IQ-^"' g cm'^ (1.48 X 
10* Mq pc~^). Recently, Kro upa et alj 1201 2f) argued that there is a den- 
sity and metallicity dependence to the IMF and proposed that the IMF 
slope above f» 1 Mq flattens to one that depends on density, qimf = 
-1.86 + O.431og(p/1O'^M0 pc~3) when p > 9.5 X 10* Mq pc-^. To 
test the observational consequences of this IMF, the suite of 1260 models run 
with the traditional Kroupa-Salpeter IMF was re-run with a top-heavy IMF 
Q-j^p = -1.86 between 1 and 100 Mq. This IMF yields e = 1.5 X 10"^. 
There were no qualitative differences in any of the results of this paper be- 
tween the two IMFs. 
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R (pc) 

Figure 1. Several properties of an example Eddingt on-lim ited star-forming 
disk calculated using the theory described in section 12.1 [ including the SFR 
(solid line), the vertical optical depth in the IR r (dotted line), the effec- 
tive temperature T^g (dot-long-dashed line), the central temperature T (short 
dashed line), the star formation surface density E» (dot-short-dashed line), 
and the gas surface density Sgas (long dashed line). 

creases the scatter of the derived relations). This selection cri- 
teria yields 132 model star-forming disks. Figure[T]shows sev- 
eral of the properties of one of these Eddington-limited star- 
forming disks. The plot clearly shows the implications of fix- 
ing Q = 1: the increasing Sgas keeps ry large enough so that 
radiation pressure dominates over nearly 200 pc. This fact, in 
turn, keeps the disk very warm with temperatures T > 100 K 
and S„ - 10^ M© yr^^ kpc^^ 

Figure [Uillustrates that the quantities that enter into the star 
formation laws (i.e., SFR, and Egas) can vary dramati- 
cally over the size of the star-forming disk. As these variations 
will be unresolved at most redshifts, we consider both the ra- 
dial average and the mass-weighted average of these variables 
when constructing the predicted KS and ES laws. For exam- 
ple, the mass-weighted average of the SFR along the disk, 
(SFR)mu,, is calculated as follows: 

^(27rrArSgas(r)//gas(r))SFR(r) Ar 

(SFR)mt„ — , 

^(2^rArEgas(r)//gas(r)) Ar 

(7) 

where Egas//gas is the total dynamical mass surface density. 
The simple radial average, (SFR), is computed as in Eq.|7]but 
with the (27rr'Ar'Egas//gas) factors omitted. 

Finally, it is important to bear in mind the limitations asso- 
ciated with this model including the lack of any dumpiness 
which is found by hydrodynamical simulations of turbulence 
supported disks (e.g., Shetty & Ostriker2012). In addition, as 
these disks are Eddington-limi ted they will may be subject t o 
Rayleigh-Taylor instabilities (iKrumholz & Thompsonll20T2h . 
and could also dr ive outflows and fee d back matter and energy 
into the galaxy (Murrav et al.l 120051: [Andrews & ThompsonI 
120111) . The effects of outflows on the structure of these disks 



is not considered here. Feedback from an accreting black 
hole is also not treated, so comparisons between the predic- 
tions and data should be limited to sources without active nu- 
clei. It is expected that the impact of these limitations will 
be mostly offset by considering radially averaged quantities, 
a wide range of model parameters, and by comparing to con- 
straints derived from observations of many dozens of individ- 
ual galaxies. 

2.2. Infrared Luminosity and SFR 

The SFR is one of the most important observational proper- 
ties of a galaxy as it gives crucial information on the galaxy's 
gas content and evolutionary state. Estimating the SFR typi- 
cally relies on utilizing one of a number of relationships be- 
tween SFR and the luminosit y of the galaxy in a specific 
wavelength or energy range ( Kennicutt & Evansl [2012). As 
most of the starlight produced in high-z star-forming galaxies 
is absorbed by dust and re-emitted in the IR, the SFR-Lttr re- 
lationship determined by Kennicutt ( 19983) and lMurphv et al.l 
(1201 II) is one of the most frequently used SFR estimators, as 
well as serving as the calibration for SFR estimates at sin- 
gle IR wavelengths (Riekeetal. 2009; Calzetti et al. 2010). 
As described bv iKennicutt (.1998b.) and .Muiphy et al., (.20111) . 
the SFR-Ltir relationship was calculated using Starburst99 
models assuming that 100% of the bolometric luminosity of a 
stellar population with a constant SFR is captured and emitted 
in the IR. This scenario is also fulfilled by the optically thick 
star-forming disk models with (ry) > 1, as ry is the vertical 
optical depth in the IR. Thus, as the SFR is known through- 
out each disk, it is interesting to determine how accurately the 
Murphy et al. -derived SFRs describe the star-forming proper- 
ties of the model disks. 

Since these disks are largely optically thick at all radii, the 
effective temperature Tcff of one side of the disk surface at 
radius r is related to half of the total flux produced by star 
formation: (JsbT^s = (l/2)eEa,c^, where a sR is the Stefan- 
Boltzmann constant (iThompson et al . 2005). Assuming that 
each annulus of the disk radiates as a blackbody and is viewed 
face-on, the IR spectral energy distribution (SED) of each 
model is calculated as 

_ 2TThc^ r°"' 2Trrdr 

where ri„ is the inner radius for each model, h is Planck's 
constant and fee is Boltzmann's constant. The Xtir for each 
model is then computed by integrating the SED from 3- 
1100 /im. This luminosity can then be compared to various 
calculations of the disk SFR, such as (SFR), SFR,„ax, and 
(SFR)„^. 

Figure |2] plots the iMurphv et"!!] (120111) calibration 
(log(SFR/M0 yr-i) = log(LTiR/erg s"!) - 43.41) as the 
solid line, SFRmax versus Ltir (open triangles), and (SFR) 
against Ltir (solid points). The SFRmax, (SFR) and Ltir 
are derived d irectly from the mo dels as described above. 
Evidently, the iMurphv et al.l (1201 Ih calormetric relation pro- 
vides a coarse approximation to the maximum SFR of each 
disk. This is a natural consequence of the assumption of 
a single SFR that underlies the original computation of the 
SFR-Ltir relationship. The maximum SFR in the disk pro- 
vides the largest contribution to the luminosity, so the esti- 
mated SFR from the Ltir relationships will be closest to the 
maximum rate. However, as shown in Fig. |2] a better (with 
~ lOx smaller scatter) SFR-Ltir relationship is found with 
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Figure 2. Radially averaged (filled circles) and maximum (open triangles) 
SFRs from Eddington-limited star-forming disks plotted as a function of the 
infrared luminosity (computed from the spectrum derived from eq. [SJ. The 
short dashed Une denotes the least-squares fit to the data while the solid line 
plots the SFR calibration published by Murphy et aU t201 h) . 

the (SFR) which is, on average, 2.2 ± 0.3 smaller than the 
values predicted by the Murphy etal. (2011) calibration at 
the same luminosity. A least squares fit to the (SFR)-iTiR 
relationship (dashed line) yields log((SFR)/M0 yr^^) = 
(0.991 ± 0.008) log(LTiR/Lo) - (10.07 ± 0.09). Thus, the 
(SFR) is used throughout the paper to quantify the SFR for 
any particular model. 

2.3. Computing the CO Line Emission 

The CO intensity for multiple rotational lines is com- 
puted for each model using the one-dimensional version 
of the molecular line radiative tr ansfer code RATRAN 
( Hogerheiide & van der TakI 120001) . I Armour & Ballantvnel 
( 20121) previouslv used RATRAN to study the CO spectral line 
energy distributions (SLEDs) of eighteen star-forming disk 
models (nine with a pc-scale burst and nine without), and a 
complete description of the calculation procedure is described 
in that paper Below is a summary of the key parameters and 
assumptions that enter into the calculation. 

The temperature and density of the star-forming disk mod- 
els are computed on a logarithmic radial grid which is input 
into RATRAN. As mentioned above the Eddington-limited 
star-forming disks are uniformly very dense, so that the gas 
is ex pected to be entirely m olecular as long as T < 1000 K 
(e.g.. lPelupessv et"ani2006h . Thus, the molecular fraction of 
the gas is set to unity for all radii, and the CO ab undance is 
set to be 10~^ of the gas number density (iKlemp erer 2006). 
Since the star-forming gas is under high pressure and typically 
has densities > 10^ cm^'', gas and dust will be well mixed, 
and the molecular kinetic temperature a nd dust tempera - 
tures are set to the central temperature T (lGoldsmithll2001h . 
Dust opacity is included in the radiative transfer calculation 
with a temperature dependent broken power-law opacity (see 
lArmour & Ballantvne..2012.) . and is scaled as appropriate for 



each model's dust-to-gas factor. The turbulent line width at 
each radius is derived from the local rms speed Cg (i.e., the 
Doppler broadening parameter is b = (2y^ii(2))^^Cs), and 
the radial velocity is racg, as defined above. Once RATRAN 
completes the radiative transf er calculation for a disk model, 
the Miriad software package (I Sault et aPl 1995 1) is used to re- 
move the appropriate dust emission from each spectral line 
and integrate over velocity to produce the velocity-integrated 
intensity Iqq in K km s^^. Using this method, CO line in- 
tensities for the J =1-0, 2-1, 3-2, 4-3, 5-4, 6-5, 7-6, 8-7 
and 9-8 transitions are predicted for each Eddington-limited 
star-forming disk model. 

3. THE STAR FORMATION LAWS 
3.1. The Kennicutt-Schmidt Relation 

Turning now to the star formation laws predicted by the 
Eddington-limited star-forming disks. Figure |3] plots the KS 
law in two ways: (E*) versus (Sgas) (left panel) and (E*) 
against (Sgas)m«; (right panel). The theoretical expectation 
is for cx Sgas/K (Eq. |6l), and, indeed, we find that the 
KS laws, calculated from these radial averages, have slopes 
close to unity and are seperated vertically by the dust-to-gas 
enhancement factor (that increases k). The scatter in the KS 
relations is smallest when plotted against (Sgas) because the 
same spatial averaging is being performed on both quantities. 
However, due to the dependence of Egas on (see eq.|4]and 
surrounding discussion), the slopes of the KS laws in the left- 
hand panel can drift away from the expectation of iV = 1. 
This effect is enhanced for larger dust-to-gas enhancement 
factors because a larger k reduces and thus Sgas does not 
decrease as fast with radius. The spatial averaging therefore 
results in a larger Sgas which then flattens the slope of the 
relations. By plotting the quanitites with different spatial av- 
eraging, the right-hand panel breaks the dependence of Egas 
on E* and shows that indeed the KS laws follow the theoreti- 
cal expectation with iV w 1. 

A random selection of galaxies will likely exhibit a wide 
range of dust-to-gas ratios that may be difficult to observa- 
tionally seperate. The dotted lines in Fig. |3]plot the KS re- 
lations found when including all the model points, regardless 
of the dust-to-gas factor The resulting slopes are — 1.2 
and 1.3 for the (Egas) and (Egas)mto panels, respectively. 
These values are very similar to the observed KS laws (e.g., 
iKennicutt & Evansl2012l) . and indicate that the common prac- 
tice of measuring the KS relation by compiling a large collec- 
tion of heterogeneous data may be hiding important clues of 
the physics of star-forming disks. Future investigations of the 
star formation laws may benefit from the analysis of a smaller 
sample of galaxies with well known properties. 

3.2. The Elmegreen-Silk Relation 

The ES law predicted from the radiation pressure domi- 
nated star-forming disk models are shown in Fig. |4] with the 
same symbols and line styles as the previous KS law figure. 
Following the current observational practice the orbital fre- 
quency is set to the inverse of the dynamical timescale at rout, 
i.e., Tdyn — ^/^{fout)- The theoretical expectation for the ES 
law is E»cx(Egasi^/K^)°'^ (eq.|5]l, which is in good agreement 
with the relations seen in Fig. |4] There is a flattening of the 
slope with the dust-to-gas enhancement factor that is reduced 
in the (Egas)mu)/Tdyn plot. The explanation for this effect 
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Figure 3. Kennicutt-Schmidt relations, defined either as (S*) vs. (Sgas) (left panel) or as (S*) vs. (Sgas)mui (right panel), predicted by the Eddington-liinited 
starburst disk models. The slope of the relations are slightly dependent on the dustiness of the gas with TV = 1.1, 0.91 and 0.84 when plotted versus (Sgas) and 
a dust-to-gas enhancement of 1 X, 5 X and 10 X , respectively. When (S,) is plotted against (Sgas)miD the slopes sue N = 1.0, 0.97 and 0.94 as the dustiness of 
the gas increases. The dotted line in each panel plots the relation when all models are considered, irrespective of the dust-to-gas factor In that case, N = 1.2 (left 
panel) and 1.3 (right panel), very similar to the KS relations measured from less intense star-forming regions at low and high redshift (e.g.. lKennicutt & EvansI 
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is the same as with the KS relations discussed above. How- 
ever, it is clear from a comparison of figs. |3]and 2] that the 
ES law has significantly more scatter than the KS relations 
(about a 3x larger rms in the (Sgas) plots; ^ 50% larger rms 
in the (Sgas) mm panels). This increased scatter is likely due 
to dividing a radial average by the orbital time at one radius, 
as well as a contribution due to variations in the rms speed 
in each model (eq. |5]l. It is notable that even when all the 
models are fit with a powerlaw, irrespective of the dust-to-gas 
factor, the measured slopes (n = 0.70 and 0.73) are still sig- 
nificantly sub-linear Thus, measurements of a sub-linear ES- 
law in high-z star-forming galaxies may provide compelling 
evidence for a high-density Eddington limited environment. 

4. THE INFRARED-CO RELATIONSHIP AND THE 
OBSERVED STAR FORMATION LAWS 

The previous section described the theoretical expectations 
for the star formation laws governing the radiation pressure 
supported disks. To test the theory, these laws must be obser- 
vationally inferred from data of a sample of galaxies, which. 



at high-z, are typically limited to photometric measurements 
such as an IR luminosity and a molecular line intensity. The 
challenge is then to turn these measured quantities into accu- 
rate estimates of and Egas for each galaxy. In this section, 
we make use of our predicted IR luminosities and Iqo for 
each disk modeQ to test how well the observationally derived 
star formation laws recover the theoretical laws presented in 
Sect. [3] 

Since an IR luminosity is correlated with the SFR 
(Sect, \2.2i and the molecular line intensity is connected in 
some manner to the total gas content of galaxy, a common 
zeroth-order measurement of the KS law is made by plotting 
the observed IR luminosity against a molecular line luminos- 
ity or intensity (Figure |5]i. The CO J = 1 — transition 
is used as the molecular line tracer in this plot because it is 
one of the most commonly measured lines; however, the den- 
sities and temperatures of these disks are large enough that 

^ We work with the velocity integrated intensity 7co instead of L'qq to 
inore closely inatch the observables expected for high-z galaxies. 
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Figure 4. Elmegreen-Silk relations, defined either as vs. (Egas)/Tdyn (left panel) or as vs. (Sgas}muj/Tdyn (right panel), predicted by the 

Eddington-limited starburst disk models. The dynamical time, r^yn, is defined as the inverse of the orbital frequency at rout, i.e., Tdyn = l/f^C^out ). As for the 
KS laws (Fig. [3), the models are divided based on the multiplicitive dust-to-gas factor with the same symbol styles as the previous figure. When (S,) is plotted 
versus (Sgas)/Tdyn (left panel) the power-law slopes are equal to n = 0.59, 0.53 and 0.49 for a dust-to-gas enhancement of 1 X , 5x and 10 X, respectively. If 
the (E*) is instead plotted against (Sgas)mu)/Tdyn (right panel) the slopes are n = 0.54, 0.54 and 0.53. The dotted line plots the least-squares fit to the all the 
points and has a slope of n = 0.70 (left panel) and 0.73 (light panel). 



all the low J CO lines are fully ther mahzed and thus have 
approximately equal luminosities (see Armour & BaUantvnel 
(I2OI2 ) for example CO SLEDs from the radiation pressure 
dominated disks). As a result, neither the slope nor the scatter 
in the Ltir-/co relationship depend significantly on the CO 
transition. 

Interestingly, the slope of the Ltir-^co relationship (1.1 ± 
0.1) is very similar to the one predicted from the theoreti- 
cal KS law (N = 1.2 or 1.3 when including all models; 
Fig. O. However, this is purely a coincidence, because in 
these warm, dense and well-mixed disks both the dust and 
molecular gas emitters are thermalized, optically thick and 
represent the same amount of disk mass. Thus, a quasi-linear 
relationship between the Xtir (indicating the dust mass) and 
- ^co (indicating the gas mass) is expected in this scenario 
jKrumholz & Thomp son 2007). We conclude that Ltir-^co 
plots are poor estimates of the star formation laws of radiation 
pressure dominated disks. 

The next step is to use the observables to estimate the physi- 



cal quantities that enter into the star formation laws, beginning 
with E,. This quantity is reasonably straightforward to calcu- 
late, as one can simply translate the Ltir into a SFR through 
one of the well known calibrations and then divide by the ob- 
served area of the disk to obtain the star formation surface 
density. Figure |6] illustrates how well that procedure recov- 
ers the actual (S,) of the model star-forming disks by plot- 
ting versu s the values obtained from Ltir, the relation- 
ship from Sect. 12.21 and the model inner and outer radii. The 
'observed' is highly correlated with (S*), and the least 
squares fit to the data (dashed line) has a slope of 1.04 ± 0.02, 
but is w 2x smaller than the one-to-one relationship. The 
dotted line is the least-squares fit to the model data if the 
Murphy et al. (2011) Xtir-SFR relationship is used and re- 
sults in a that is approximately correct at low (S*), but 
becomes progressively too small at large (S,). These offsets 
from the actual (S*) is a result of approximating the aver- 
age with ((SFR))/area or (SFRmax/area), but are relatively 
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Figure 5. The Ltir and Iqo relationship obtained from the Eddington- 
limited star-forming disk models. The dashed line shows the least-squares 
fit to the model data with a slope of 1.1 ± 0.1. If this plot is interpreted 
as an observational analogue of the KS law, then the slope is in reasonable 
agreement with the theoretical expectation of 1.2 or 1.3 (Fig.[3). However, 
this quasi-linear relation is simply a result of the molecular Une and dust 
emission tracing the same gas population, and is not an accurate estimate of 
the underlying KS law. 
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Figure 6. A comparison between the (E*) predicted by the Eddington- 
limited starburst disks and the one calculated from the 'observed' Ltir, a 
Ltir -SFR correlation, and the inner and outer radii of the disk. The solid 
line plots the one-to-one correspondence between the two quantities, the dot- 
ted line is the relationship found when the lMurphv et al.l 8201 1) LxiR"SFR 
cahbration is used to derive (S*), while the data points an d the dashed line 
result from the Ltir-(SFR) correlation described in Sect. 12.21 This figure 
indicates that, due to the spatial variation of the stai' fonnation rate, any sim- 
ple estimate of (S» ) will be inaccurate by a factor of < 2. 



minor. Most importantly, there is a good linear correlation 
between and the observational estimate, which indicates 
that the translation from Ltir results in a reasonable estimate 
of (E,). 

As described in Sect. [T]the conversion factor aco between 
Iqo and Egas is not well determined except for a small num- 
ber of galactic star-f orming region s and s ome local ULIRGs 
(e.g., Downes & SolomonI 119981 : ISolomon & Vanden Boul 
120051: Naravan anetal.1 120121) . In our case, since we have 
both Iqo and (I]gas)mM> for every model, we can calculate an 
unique value of aco for each star-forming disk (calculating 
aco with (Egas) instead of (Egas)m«i results in more scatter, 
so we focus on those computed with the mass-weighted av- 
erage). Figure |7] plots the relationship between the resulting 
aco ™d Iqo and indeed finds a reasonable anti-correlation 
(R — —0.64, where R is the linear correlation coefficient) 
with a slope of —0.56 ± 0.06. The median value of aco is 
3.56. As the CO intensity is correlated with the SFR (Fig.|5), 
this plot also indicates that aco is anti-correlated with the 
SFR. The slope of the aco-^ no anti-correlation is s teeper 
than the one of —0.32 found bv lNaravanan et aLl (120121) from 
hydrodynamical simulations of galaxies. Similarly, the am- 
plitude of the do tted line in Fig. |7]is > 2x larger than the 
[Narayanan et al.l(l2012l) relation. As the physical scenario and 
calculation of the CO intensity performed by Narayanan et al. 
(|2012) differs significantly from our setup the similarity in 
the results is striking. The relatively high gas temperatures 
found in the Eddingt on-Umited star forming disks will tend 
to dri ve aco lower dMagnelli et alj l20T2t iNaravana n etaTI 
I20T2I) . which may explain the steeper slope in Fig. |7] The 



larger amplitude may be caused by the lack of dumpiness 
in our star-forming medium which, when combined with the 
optically-thick nature of the CO emission, would result in 
the CO line intensity being a worse tracer of the total gas 
mass than in a clumpy medium where the CO intensity is 
directly related to a larg e number of smaller clouds. How- 
ever, ^apadogoulosetan I2OI 2) recently argued that the cen- 
ters of ULIRGs (which may be radiation pressure dominated: 
[Thompson et al.ll2005h are unlikely to be comprised of indi- 
vidual molecular clouds and that the small values of aco in- 
ferred from these environments may be erroneously low, in 
agreement with the results presented in Fig. Q 

It is apparent from Fig. Qthat there is no single value of 
aco for a given Iqo ; indeed, the scatter in aco at one Iqo 
can reach a factor of ten. Using the best-fit aco'-'^CO rela- 
tion will than propagate this large scatter through to the esti- 
mated Egas (as seen in Fig. [8]). The values of Egas derived 
from aco are only weakly correlated (i? = 0.55) with the 
(Egas)m«) calculated directly from the models. This poor re- 
lationship between Iqo and Egas will substantially compli- 
cate any attempt to recover the star formation laws from the 
observational quantities. 

Finally, Figure |9]presents the KS and ES laws of the radia- 
tion pressure supported star-forming disks derived solely from 
the observational quantities Ltir and Iqo ■ To guide the anal- 
ysis and comparison with the theoretical expectations (figs. [3] 
andlH, the models are separated by their dust-to-gas enhance- 
ment factor, although this may be impractical for a real ob- 
servational sample. Recall that the Egas computed from Iqo 
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is an estimate of (Sgas)miu, so Fig. |9]should be compared to 
the right-hand panels of Figs. [3] and |4] As expected, the scat- 
ter of the 'observed' KS laws (rms« 0.2) is larger than the 
ones found from the theoretical relations (rmsw 0.1). More 
striking are the differences in the correlation coefficient and 
slopes. The observationally derived KS laws have slopes of 
= 1.3 ± 0.3, 0.8 ± 0.2, and 0.5 ± 0.2 and coiTelation 
coefficients of i? = 0.6, 0.5 and 0.4 as the dust-to-gas en- 
hancement factor increases. In contrast the theoretically ex- 
pected slopes are 0.94-1.0 and have correlation coefficients 
of > 0.9. Collecting all the models together yields a KS slope 
of iV = 1.4 ± 0.2 (with R = 0.5), surprisingly close to the 
canonical value, but, of course, yielding an entirely mislead- 
ing view of the underlying physics. Thus, the KS law gov- 
erning the Eddington-limited star-forming disks is nearly un- 
recognizable due to the very poor conversion between Ico 
and Egas- Our data indicate that the use of a higher den- 
sity CO transition (such as CO J = 3 — 2) would not im- 
prove the conversion as this line is thermalized in these very 
warm disks. A high- J transition from a large density trac- 
ing molecule, such as HCN, may provide the best probe of 
Eddington-limited disks (although this needs to be confirmed 
by subsequent computations). 

In contrast to the KS law, the 'observational' ES laws 
(right panel of Fig. |9]l have similar correlation coefficients 
(R w 0.8) and about the same amount of scatter (rmsw 0.2) 
as the theoretically expected relations (right panel of Fig. |4|i. 
However, the slopes of the different ES laws are all n > 0.7, 
inconsistent with the expected slope of n — 0.5. The slope 
for the entire ensemble of models is n = 1.2 ± 0.1, which, 
much like the KS law, leads to an erroneous description of the 
correct ES law. The steep slopes are a result of the uncertain 
Q^CO-^CO relation. If, for example, a constant aco — 3.56, 
was used to estimate Sgas, then the slopes flatten to n = 0.4- 
0.7 for the individual dust-to-gas factors, and n — 0.8 ± 0.1 
when including all the models. The se re sults are consistent 
with the theoretical ES law from Sect. 13.21 Applying this con- 
stant aco to the observed KS law results in a slope that is 
much too flat (N = 0.2-0.6). The slopes of the observed ES 
and KS relations are therefore very sensitive to the assump- 
tions placed on aco; however, the low rms and large corre- 
lation coefficients found in the observed ES law are robust to 
the assumptions on aco^ indicating that, out of these two rela- 
tionships, the ES law may provide the most promising means 
to investigate SF laws at high redshift. 

5. DISCUSSION & CONCLUSIONS 

The purpose of this paper is to, one, elucidate the KS and 
ES laws of a Eddington-limited star-forming disk (including 
the impact of spatial averaging), and, two, to determine how 
well the theoretically expected star formation laws can be 
recovered by the standard conversions from the observables 
(Ljfj and Ico)- To explore the effects of spatial averaging 
and how these disks might appear to telescopes, 1260 star- 
forming disk models, spanning a wide range of size and mass 
scales, were run to produce a bank of test data. From this 
model suite, 132 of the largest and most optically thick disks 
were selected to best describe the properties of Eddington- 
limited starbursts that may be embedded in the centers of the 
gas-rich star-forming galaxies. 

The analytical model of a radiation pressure suppo rted star- 
forming disk developed by [Thompson et al.l ( l2005h predicts 
a KS law with w 1, and an ES law with n w 0.5, both 
significantly flatter than the canonical values {N sa 1.5 and 
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Figure 9. The KS (left panel) and ES (right-panel) laws of radiation pressure supported star-forming disks derived from the Ltir-^CO relationship (Fig. |5) 
as described in the text. As in the plots of the theoretically expected star formation laws (Figs. |3]and|4), the models have been separated by the value of the 
dust-to-gas enhancement factor The slope of the relations are dependent on the dustiness of the gas with the KS relations yielding N = 1.3, 0.8 and 0.5 for a 
dust-to-gas enhancement of 1 X, 5x and 10 X, respectively. For the ES relations the slopes are n = 1.1, 0.8 and 0.7 as the dustiness of the gas increases. The 
dotted line in each panel plots the relation when all models are considered, irrespective of the dust-to-gas factor In that case, TV = 1.4 ± 0.2 (KS law) and 
1.2 ± 0.1 (ES law). Overall, the laws derived from observations are poor representations of the actual star formation relations that govern the disk. 



n « 1). An important aspect of these star formation laws is 
the dependence on the opacity, or, equivalently, the dustiness 
of the star-forming gas. Indeed, the star formation process in 
a ra diation supported disk strongly depends on the dust opac- 
ity ([Thompson et al.i r2005), and, when spatially averaged, can 
even impact the slopes of the observed star formation laws 
(Figs. [3]and|4]i. Another important effect of the gas opacity 
is the vertical separation of the star formation laws based on 
the dust-to-gas enhancement factor (Figs.|3]and|4|i. A hetero- 
geneous sample of Eddington-limited star-forming disks may 
then appear to have a significantly steeper slope than what is 
governing the disks and would lead to the wrong conclusion 
about the underlying star formation relations. 

To study how these disks might appear observationally, 
molecular line radiative transfer calculations were performed 
and CO line velocity-integrated intensities, /co. were de- 
rived for all 132 star-forming disk models. Combining these 
results with the infrared luminosity, Ltir, calculated from 
the SED of each model, allowed a quantitative test of de- 
riving the star formation laws from observational data. The 



calculated Ltir and Ico are approximately linearly corre- 
lated with each other (albeit with a substantial scatter of 
rms= 0.4; Fig. |5]l, as might have been expected from the 
theoretical KS law. However, in the Eddington-limited star- 
forming disks, the gas density and temperature are always 
above the critical density and excitation temperature of the 
first several CO rotational transitions. These lines are there- 
fore completely thermalized and are tracing the same gas as 
the infrared emission, resulti ng in a near-linear relationship 
(e.g., iKrumholz & ThompsonlEOOTt [Narayanan etaH 1^008). 
Indeed, the densities and temperatures of these disks are 
so large that a linear relationship is also expected for the 
low J transitions of molecules with higher critical densities 
(such as HCN or HCO''') , as seen in some ULlRGs (e.g., 
iGarcia-Burillo et aDl2012l) . While such a linear relationship 
may indicate that the star-forming region could be radiation- 
pressure dominated, i t does not give any informa tion on the 
underlying KS law (Kru mholz & T hompson '2007). Obser- 
vations of high- J transitions from these molecules may be 
needed to determine the KS and ES laws from the observed 
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IR luminosity and molecular line intensity. 

If the ixiR-Zco relation cannot be trusted as a measure- 
ment of the star formation laws of Eddington-limited disks, 
then these quantities must be translated to and Sgas- Fig- 
ure |6] shows that the use of a Ltir-SFR conversion and the 
outer radius of the disk results in a (S*) that is within a fac- 
tor of 2 of the correct value. On the other hand, aco^ the 
phenomenological conversion factor between Iqo ™d Sgas, 
exhibits significant scatter (about an order of magnitude for a 
given /co; Fig- 13 and yields a very weak correlation between 
the estimated Egas and the true (Egas)mtu (Fig.©. This poor 
conversion from Iqo to Egas causes the resulting KS laws 
(Fig.|9](left)) to have only a weak correlation between E, and 
Egas. The slopes of this KS relation vary strongly with the 
dustiness of the gas, reaching as low as = 0.5 ± 0.2 for a 
dust-to-gas factor of 10. The slope of iV = 1.4 ± 0.2 found 
from including all the models is consistent with the canonical 
slope, which perfectly illustrates how the observed KS law 
can distort and mislead the interpretation of star formation 
laws. 

The situation is improved somewhat when considering the 
ES law, calculated by dividing the /co -derived Egas by the 
orbital time at rout (Fig. |9] (right)). In this case, the observed 
ES laws retain the large correlation coefficients expected from 
the theoretical predictions, but, as with the KS law, the slopes 
are significantly steeper (n > 0.7) than what the theory pre- 
dicts (n « 0.5). Similar to the KS law, the slope found from 
the ensemble of models {n — 1.2 ± 0.1) is close to the canon- 
ical linear slope, leading to yet another erroneous conclusion 
about the underlying star formation law. 

These experiments all point to the same problem: the aco 
conversion factor is too blunt and inaccurate an instrument 
to study the star formation laws in radiation pressure domi- 
nated disks. In many ways this is unsurprising, as the aco 
factor has only been shown to accurately trace Egas in lo- 
cal Galactic molecular clouds up to Egas~ 100 Mq pc~^ 
fKennicutt & Evans! 120 12) . The Eddington-limited disks are 
at such high surface densities that the medium may not be 
accurately described as a series of molecular clouds, but 
rather a continuous structure with a near unity molecular frac- 
tion. In their study of the CO SLEDs from these disks, 
lArmour & Ballantynd (|2012) noted that in many ways the 
thermalized molecular emission from these disks acts like a 
stellar spectrum, with a SLED shape that depends on exci- 
tation temperature and a normalization that depends on the 
internal energy generation, i.e., star formation. Therefore, a 
large, massive disk may produce a similar Iqo as a smaller, 
less massive disk that has a larger (SFR) (due to, e.g., a larger 
opacity). The CO emission is largely immune to to the to- 
tal gas mass in each disk, and therefore the aco is a poor 
tracer of Egas. We do find that aco is invers ely related to Igo 
and ( SFR) , in agreement with other authors ( Narayanan et al.l 
120121) . but with an amplitude ss 2 x larger This may be an 
effect of the high density and uniform molecular medium of 
these disks (iPapadopoulos et al. 2012) . but this result should 
be checked with more reaUstic two- or three-dimensional sim- 
ulations. 

Given the difficulties in using aco^ how can we make 
progress in studying star formation laws, especially when 
Sgas> 10"* M© pc"2? Clearly the traditional method of CO 
rotational lines will not be particularly useful in this regime, 
so observations should focus on molecular species with higher 
critical densities, such as HCN. Recent results using HCN ob- 



servations of ULlRGs that reach into this high- Egas range do 
give KS and ES slopes consistent with the theoretical expec- 
tation for Eddington-Umited disks, suggesting that radiation 
pressure dominated regions may b e playing a role in these 
sources toarcfa-Burillo et ani2012h . Unfortunately, the con- 
version factor of HCN and other higher density tracers is as 
uncertain as aco- Observational campaigns focusing on mea- 
suring these conversion factors over as wide a range of Egas as 
possible will be necessary in order to confidently apply them 
for ULIRGS and at high redshift where large gas surface den- 
sities will be more common. Interestingly, the ES law ap- 
pears to be more robust to changes in opacity and maintains 
a strong correlation even with the use of aco (although the 
slope is still incorrect). This result implies that, at least for 
the Eddington-limited star-forming disks studied here, mea- 
surements of the ES law will be more useful in determining 
the underlying star formation relations than the KS law. Fi- 
nally, the star formation laws of radiation pressure dominated 
disks are strongly dependent on the gas opacity, causing a ver- 
tical separation in the KS and ES plots. Fitting a line through 
the models without controlling for the differences in the dust- 
to-gas ratio results in a much steeper slope than the under- 
lying star formation law. Thus, future observational work at 
these values of Egas would benefit from analyzing samples 
of galaxies with metallicity or dust-to-gas ratio estimates in 
order to control for this potentially misleading effect. 

In summary, the KS and ES star formation laws for radia- 
tion pressure supported star-forming disks have flatter slopes 
(N w 1 and n w 0.5) than what is expected for star formation 
at Egas< 10'' Mq pc^^. These slopes are relatively robust to 
spatial averaging over the disks, but the strong dependence on 
the opacity causes the relations to be vertically offset based on 
the dust-to-gas ratio. Thus, an erroneously large slope for the 
KS and ES laws could result if the opacity dependence is not 
recognized. Attempting to recover these star formation laws 
from the predicted Ltir and Iqo failed because of the poor 
translation from Iqo to Egas at these high densities and tem- 
peratures. The 'observed' KS and ES laws have slopes that 
differ significantly from the theoretically expected relations, 
even when controlling for the gas opacity. Progress in iden- 
tifying and studying the star formation laws at these values 
of Egas can be made by moving to well calibrated high crit- 
ical density molecular tracer, selecting a sample of galaxies 
with metallicity estimates, and focusing on the ES relation, as 
its slope seems less affected by variations in metallicity and 
spatial averaging. 

The authors thank D. Narayanan for comments on a draft 
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